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We outline the nonlinear transformation in the path integral representation for 
partition function of O(N) symmetric oscillator systems bringing their duality to 
certain one-dimensional oscillators with unstable potential shapes. This duality 
transformation realizes the equivalence between spectra of a Hermitian and a non- 
Hermitian Hamiltonians. It is exploited to calculate exactly the spectra of quantum 
rotators . The zero-temperature limit is considered and the dual representation of 
n-point correlators of coordinate operators is obtained. 

1 Introduction 

The method of equivalent oscillators (EO-transformation) elaborated in [T] allows one to 
find the correspondence between quantum oscillator systems with attractive and repulsive 
anharmonic interactions. This equivalence can be realized when the coordinate paths for 
the repulsive interaction are taken along a curve in the complex space jH [2] . Recently 
this kind of relationship between Hermitian and non-Hermitian oscillator systems has 
attracted much attention 0, H] in the search for examples of PT-symmetric non-Hermitian 
quantum Hamiltonians with real spectra [2j [5j EJ [7] . 

In our paper we present the derivation of EO-transformation using the phase-space- 
like path integral, find a novel, general representation for dual QM system and apply it 
to the partition function of anharmonic oscillators and of quantum pendulum in order to 
calculate their spectrum by path-integral means. The n-point coordinate correlators are 
examined and their dual analogs are obtained. 

The basic EO-transformation evaluates the partition function Z(/3), f3 = 1/kT for a 
iV-dimensional harmonic oscillator with variable complex frequency Cj 2 {t) = u 2 + io~(r) 
assumed to be periodic, cu 2 (0) = cu 2 (/3) , 
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where (q) = (q±, . . . , q^r). From the structure of the action in ([I]) it is seen that the 
imaginary part of the frequency may serve to realize the Legendre transformation to the 
description dual to the O(N) invariant coordinate q 2 . The real part of the frequency 
(I> 2 (t) could also depend on Euclidean time but in this case we assume its positivity 
W, Re u) 2 {t) = uj 2 + u(t) > 0. Let's skip the variable part u(t) for a while as it can be 
always embedded by shifting a — > a — iv. Formally, after path integration one comes to 
the functional determinant of the differential operator — d 2 + uj 2 , 
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where the constant C is chosen so that for the constant u = uj the partition function is 
saturated by the discrete spectrum of the harmonic oscillator, 
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This analytic expression is not available when the function <j(t) is variable in time Instead 
it is given by the functional nonlocal in cx(r), 
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The temperature Green function G^t) is defined to fulfill the periodicity in Euclidean 
time G P (Q) = G (P), 
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Our task is to perform the local change of the variable <t(t) so that to find explicitly this 
functional in terms of new variables. 

Let's introduce the nonlinear transformation, 



Co 2 (t) = uj 2 + w(t) 



-u[t) + u (r), 
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with u(t) being periodic, u(0) = u(j3) . One can show pQ that in order to provide the 
one-to-one correspondence between a(r) and u{r) one has to fix the sign of Re u(r), for 
a definiteness, we adopt Re u{t) > 0. 

The original action can be re-expressed as 
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after exploiting the periodicity of involved functions. 

Consequently the calculation of the functional determinant (J4]) is equivalent to evalu- 
ation of the Jacobian of the following change of (periodic) variables, 

'q[r) + u(T)q(r) = &(r) + wft(r); (8) 
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where the periodic Green function for inverse first-order differential operator is given by 

Gf{j-T X ) = expt-^r-rOl^r-n) 
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with 9{x) being a step function. Eventually in these variables the Jacobian reads, 
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and therefore the partition function becomes a local function of / cItu(t) 
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This is the master formula of the EO-transformation. We remark that in the final ex- 
pression the parameter u is hidden in the definition of the variable u(t), namely Re 
(—it + it 2 ) = u 2 , Im (— u + u 2 ) is arbitrary, according to Eq. ([6]) . The overall sign of 
the denominator in (fTT!) conforms to the condition Re u > . For the opposite choice Re 
u < the sign of the integral in (llip must be changed u — > — m . One can use Eq. (TTTT) 
directly to find the partition function for a temperature dependent oscillator frequency, 
just fitting u{t) instead of <t(t). 

However we intend to use the functional of er(r) to perform the Legendre transforma- 
tion from the multidimensional O(N) symmetric QM to a dual one-dimensional QM with 
the same energy spectrum. Namely let us consider the spherically symmetric oscillator 
potential uj 2 q 2 + V an (q 2 ) with an anharmonicity V an bounded from below and substitute 
the quadratic invariant <f 2 (t) by an independent variable j(t) . In the path integral for 
partition function it can be realized by means of the auxiliary Lagrange multiplier cr(r), 
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Certainly the usage of the function cr(r) as a new dynamic variable in the path integral ([I]) 
leads to a non-local Quantum Mechanics of this "observable" due to (jSJ) and (J4]). Mean- 
while the EO-transformation ([6]) makes this QM local in terms of the complex variable 
u{t) as is encoded in Eq. (TTTT) . 

To complete the EO-transformation it remains to calculate its Jacobian in the integral 

(H2D, 
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This Jacobian for periodic boundary conditions has the same structure as the Jacobian 
in the change of coordinates (|S]) and can be calculated in a similar way, 
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The infinite constant C must be absorbed by the redefinition of the path integral normal- 
ization. The overall sign in (1141) again conforms to the condition Re u > . 

We remark that in the limit u — > this Jacobian vanishes as the relation (I13p between 
5a(r) and 5u(t) in this limit is not invertible. Instead for Re u > it is invertible and 
the Jacobian (fill does not vanish. 



2 Dual potential systems 

Thus the entire duality transformation relates the N-dimensional O(N) symmetric system 
with a potential oo 2 q 2 + V an (q 2 ) to the one-dimensional quantum system with complex 
momenta, 

Z p = J V N q(r) 
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where the complex momentum u runs along the trajectories with Re (— w + u 2 ) = u 2 , 
and Re u > 0. We notice that although the variable j(t) replaces the positive invariant 
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q 2 (t) it can be also treated as an unconstrained real variable running along the entire 
axis according to lfP2l) . To make it consistent one can define V an = V^ n (|7(r)|). 

This representation contains the canonically conjugated variables 7(1") and w(r) with 
the kinetic term 7(t)u 2 (t) mixing both conjugated variables. Therefore the path integral 
suffers from the ordering problem [H [9] . As it was pointed out in [1] the Weyl ordering 
suits well to justify the periodicity conditions and the vanishing of boundary contributions. 

Expanding the function, 

<S>(z)=Y,a<i N) exp (-(- + /-!),); 
1=0 z 

4»)^ (JV+2 i_ 2) (yl±i__| forJV>1 . (16) 

= exp(-z) + exp(— -z) for N = 1, 

i.e. a^p = di = 1; a\ i\i >2 — 0, (17) 

one obtains the statistical weights (degeneracy numbers) in the decomposition of the 
partition function into the sectors with a given angular momentum I, 

OO 
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x exp | — y dr — u(t)j(t) — 7(r)n 2 (r) 
o 

W 7 (r) + Kn(7(r)) + (y + i - l)«(r)] }. 

The weights correctly reproduce the degeneracy numbers for N > 1, for instance, in 
the case of 0(2) symmetry a\ = 2, for 0(3) symmetry a\ \ = 21 + 1 etc. 

The further evaluation may be two-fold following a choice of integration domain for the 
variable j(t). According to the ansatz ( fl2|) one can use either constrained, non-negative 
7 = p 2 > or unconstrained — oo < 7 < 00. For the first choice one may safely deform 
the complex path u(r) = p(r)u(r) holding Re u(r) > and therefore not intersecting 
(P \ 

the poles of $ / dru(r) in the path integral ffl5l) . The further shift u(r) — ► u(r) — p 
\o / 

leaves the functional $ invariant due to periodicity p(0) = p(f3). Finally, for N > 2 one 
can complete this evaluation by contour deformation 

-/ n ~ / \ N + 21-2 
u(r) - u(r) + 2p(r) 
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which again does not cross the poles of $(2) just increasing the real part of complex 
momentum variable. One obtains, 
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(N + 21 -2] 2 -l,, 
WV (t) + V + 1 16p2(r) ; ] }, (19) 

which is nothing but a partition function for the reduced radial Schrodinger equation 
for ^-dimensional particle in the spherically symmetric potential for a given angular 
momentum /. The extra contribution — l/16p 2 in the centrifugal term of (IT^|) is due to 
correct resolution of ordering problem [I] in the chain of above transformations . The 
representation ffl9|) is justifiable for N > 2 . However it is not valid for N < 2 due to 
the fact that the EO-transformation is unambiguously defined only for Re u > and for 
N < 2 the shift of u(t) required to calculate the Gaussian path integral crosses the line 
Re u = confronting the pole problem in <&(z) of Eq. ( JT3T) . 

For the second choice of the unconstrained coordinate variable j(t) , in a special case 
of quartic anharmonicity V an (q 2 ) = X(q 2 ) 2 , one can instead integrate in 7(1-) to derive, 

{u{r)f 



S[u, 7; A] => S[u; A] = J dr 
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This representation is not subject to any bounds on iV and is valid also for N < 2 . 

In perturbation theory one expands [1] around a zero-order stationary configuration 
along the imaginary axis u{r) = u + i2y/Xv(r) with a real fluctuation v{r) . The latter is 
achieved by the deformation of complex paths u(r) , 

P 

S[v; A] = / dr\(v(r)) 2 - 4Aw 4 (r) + 8iVA^ 3 (r) 

(21) 
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+4u 2 v 2 (r) + (- + I - l) (u + i2V\v(r)) 



As a result, within the perturbation theory, the set of anharmonic oscillators with the 
complex quartic potential unbounded from below is discovered to describe the same energy 
spectrum as the original real O(N) symmetric Hamiltonian in ffl5l) with the potential 
bounded from below. 

This equivalence was first found for the dim-2 oscillator [TUl EEE] in the singlet sector 
I = in several (up to 15) orders of perturbation theory (see the discussion in [12j) and 
proved for arbitrary I in [1] . It was confirmed later in [13] using the technique of differential 
equations and integral transforms. For one-dimensional oscillator N = 1 with u = a 
similar equivalence has been derived recently in [21 H] . In the latter case evidently two 
branches of trajectories u with different signs of u coalesce. The case N = 1 needs a special 
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comment because one finds in (fT8|) two contributions into the partition function with the 
Lagrangians related by reflection u <-> —u [14J. One comes to its right interpretation in 
the perturbation theory, say, for the anharmonic oscillator (|2"T|) . Indeed, the lowest order 
in the anharmonicity A is saturated by two harmonic oscillators with double frequencies 
2uo and energy shifts +u/2 . Therefore two spectra u(2n + 1 ± 1/2), n = 0, 1, . . . , jointly 
form the entire ladder of levels of the initial oscillator u(n + 1/2). 



3 Duality for quantum pendulums 



Let us now use the EO-transformation for the analysis of a-models which in QM corre- 
spond to quantum pendulums. The simplest one is defined by the constraint (|12p for a 
fixed j(t) = R 2 {t) embedded in the path integral for partition function (pQ) . It repre- 
sents a quantum pendulum with variable radius R(t) . The EO-transformation eventually 
produces the decomposition analogous to ([TBI) and ffTSj) . 
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where the last line is obtained for a conventional rigid rotator with ultralocal quadratic 
actions admitting the exact calculation of the energy spectrum E t = (N + 21- 2) 2 /16i2g 
after the contour deformation u — > iu. This spectrum dependence of I has been first found 
in [15] for N = 3 with a different value of the ground state energy. The latter discrepancy 
is explained by a different choice of pendulum coordinate variables and correspondingly 
by a different choice of ordering the momenta and coordinates. Our choice is correctly 
normalized to get a zero ground state energy for N = 2 when the quantum pendulum 
is characterized by the Schrodinger Hamiltonian —R 2 d\ in the polar-angle coordinate 
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< <p < 2tt. In this case E\ = 1 2 /4Rq. In the general case of rotator with temperature 
(Euclidean time) dependent radius R m in < R{ T ) < R-max one could normalize 
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L o R(t) 



dr In 
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4 Dual form of correlators 

The partition function contains all the information on the energy spectrum and not much 
knowledge of wave functions. Therefore it makes sense also to analyze its response to 
external classical sources ff(r) with modification 



S(q,ff) = S(q,0) -J drq ■ rj (r) 



(23) 



Conventionally the expansion of the partition function ffl5|) with the action ff23l in pow- 
ers of external source f]j{r) leads to the correlators of coordinate operators at different 
Euclidean time, 
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The main ingredient of this generating functional is given by modification of Eq.([T]), 
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where we recall that Co 2 = uj 2 + ia . When applying the EO-transformation er(r) — > u{t) 
one comes to the following form of the Green function, 



Gf(r 1 ,r 2 \a) = (r\ 
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The manifest form of this Green function is easily obtained with the help of convolution 
of two first-order Green functions of (Q type, 
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Once the full information on energy spectrum and wave functions of the Hamiltonian 
can be derived from correlators it is sufficient to analyze them at zero temperature, i.e. 
at (3 — > oo . For this limit let's symmetrize the integration interval < r < (3 =>• 
—(3/2 < r < +13/2 and u{r — — > u(t) to get the entire axis in the zero-temperature 
limit. Certainly in this limit only lowest, ground state energy level contributes and the 
generating functional may be well provided by path integral over coordinate functions 
vanishing at the infinity, q(r) — > when r — > ±oo . Obviously Eqs. f|T5|) and f|T9|) are 
saturated with leading contribution for / = as Re u > 0. The formula for correlators 
( |25l) also survives with the limiting Green function, 

oo 

Grg(7i,T 2 |«)= / dreiT-nWT-Ti) 

— oo 

oo 

Xexp(- J dT'u{T')J{T',T,T X ,T 2 )), (28) 

— oo 

J( 2 V, t, n, r 2 ) = 0(r - r')[^(r' - n ) + 0(r' - r 2 )]. 

Thus the two-point correlator (propagator) of coordinate operators is given by averaging 
a "holonomy" exponential along the contours connecting two points. These contours are 
to provide the decreasing of exponentials as Re m > 0. 

After passing through the EO-transformation technology one comes to the generating 
functional of correlators for an anharmonic potential, 



Z( V ,v) = J 
x exp | - 



Vu(t)Vj{t) 
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u(r)j(r) — 7(t)m 2 (t) 
'N 



+ (^ 2 + u(t)Ht) + V an (j(r)) + (-- l)u(r) 

^ °° °° N 
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— -( 

Z(0) 



-oo — oo 

1, 
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where we have introduced also the real part of harmonic frequency uj 2 — > uj 2 + z/(r) as an 
external source which allow to generate independently the correlators of singlet operators 
(q) 2 (see the beginning of this article). This option will produce a number of identities 
between local and non-local description of the singlet sector (see below). 

From (1291) the n-point correlators are described with the help of a specific external 
source in the Euclidean Lagrangian, 



■ ■ ■ qj(2n)[T2n y 



1 " 

on n E/ 

1=1 ft}U{m,} 



x J rfr/(5 j(fe;)j(mi) 0(r/ - T kl )9(rl - r m; )H(r, r'), 

— oo 

Vu{t)V^{t) 



S({r fc },{r/}) 



const 



x exp 



{ 



dr 



(30) 



u{t)^{t) — 7(r)-u 2 (r) + uj 2 ^{t) 

+V an {i{r)) + l)u(r) + J^(t)u(t) , 

where dim{^} = dim{m;} = n, 

{ki} U {mi} = Perm{l, • • ■ , 2n}, {h} n {mi} = 0, 
and the external current includes all contours linking triples of points r/ with and r m% , 

n 

J {2n \r) = £0(r/ " t)[0{t - r kl ) + 9{r - r m )\. (31) 



i=i 



We remark that at any point r this current takes integer values. In the radial coordinate 
representation for S({rfc}, {r/}) analogous to (1T91) one obtains the action, 



5(p|{r fc },{r/}) 



^(r)+wV(r) 



+V(p 2 (r)) + 



(A^-2 + 2j( 2 ")(r)) 2 - 1 



16p 2 (r) 



(32) 



Thus the current j( 2n )(r) replaces, in fact, an eigenvalue of angular momentum depending 
on a configuration of piece-wise contours linking r{ with and r m; . 

For the case of quantum pendulum with constant radius (1221) one finds the kernel of 
n-point correlator in the following form, 



(33) 



exp 



dr 



J^ 2n \r)(N -2) + (J (2n) (r))' 
4R 2 
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which provides the n-point correlator (!30|) with an integral representation corresponding 
to generalized hypergeometric functions. Thus according to the common understanding of 
what is integrable model, O(N) symmetric quantum pendulums are completely integrable 
in dual variables. We find it promising in establishing ODE/IM relationship between one- 
dimensional differential equations and low-dimensional integrable models 
One can easily find a similar action for the anharmonic oscillator fl20|) . 

(u(r)f {u\r) - ^ - u{r)f 



S[u; A] = J dr 



4A 4A 



+( 



N +J^(t)-1 + ^)u(t)1 (34) 



2 2A 
supplemented with the external source u(t) . With its help one can find further non- 
trivial relations between " loop" and local representations of correlators for coordinate 
operators. For instance, the v.e.v. of the singlet operator ^(<f) 2 (Ti)^> can be obtained in 
two ways: first by variation in v(t), 



r]=v=0 



-L(u(t)-u 2 (t)+u; 2 + v(t)) u , (35) 



2A 

second, as a two-point correlator at coinciding arguments T\ = r 2 



»T 00 

\ 2 



5 2 ff(Ti) ' " ' ri=v=Q 

J dr'exp -2 / dru{r) (36) 



V 



Thereby the nonlocal "loop" integral (|36j) can be identified with a local differential func- 
tional (1351). 

It is not difficult to establish a similar relation for quantum pendulums checking that 
(q)\ri))=R 2 . 



Let's make few concluding remarks: 
We have proved the efficiency of the EO-transformation not only in making different 
dynamical systems equivalent in their spectra but also in providing the dual form of 
correlators which guarantees the relationship between wave functions. The correlators for 
dual oscillators have been recently discussed in [3] and one can compare our results to see 
the progress in their description . 

We notice also that an analogous method is applicable to transform the evolution 
operator for the harmonic oscillator with a time dependent frequency. 
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